Computer simulation of structural properties of dilute aqueous solutions of argon at 

supercritical conditions. 
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Abstract 

Computer simulation studies of aqueous solutions of argon are performed from ambient to supercritical conditions by using a 
recent polarizable potential model (BSV) and the non polarizable simple point charge extended (SPC/E) model. At T = 673 K 
we compare the water-solute pair correlation functions of the argon-water mixture with recent experimental results obtained 
from neutron scattering experiments. The comparison shows that the introduction of the polarizable effects decreases the 
solute-water repulsion and improves the agreement with the experiment at supercritical conditions. In particular we find that 
the water-solute structure predicted by the polarizable model is in good agreement with the experiment. 



PACS numbers: 61.20.Ja, 61.20.-p, 61.25.-f 



I. INTRODUCTION 

The problem of solubility of non polar solutes in water 
has become of large interest for the applications in many 
chemical processes in industrial technology and biochem- 
istry.Era At ambient conditions the solubility of non polar 
gases in water is some order of magnitude lower than 
in other solvent like liquid hydrocarbons, for this reason 
they are usually considered as hydrophobic solutes. The 
solubility however increases steeply above the water crit- 
ical point, T c = 647.13 K and P c = 220.55 bar. Most 
applications involve the use of water at supercritical con- 
ditions, where it, becomes a good solvent for several non 
polar materials.^ 

In fact the curve of solubility of non polar gases at 
increasing temperature shows a very peculiar behaviour, 
starting from room temperature it decreases and reaches 
a minimum, then it sharply increases toward water su- 
percritical temperature, .where complete miscibility is 
reached for some system.Q It is also been observed that 
the solwbility of rare gases is larger for larger size of the 
soluteBtl The behaviour is determined by the interplay 
between energetic and entropic contributions to the free 
energy of solvation. According, to the simulation work 
done by Guillot and GuissaniEuwith the simple point 
charge extended (SPC/E) modelQ for water the energetic 
term favors the solubility of larger solutes while the en- 
tropic term depresses the solubility at increasing size. 

The size of the solute seems also to be connected to the 
different behaviour of the mixtures of water and iWjblc 
gases at high pressure at the gas-gas critical point .Em 

Structural and computer simulatioiiElEj studies on 
aqueous solutions of rare gases indicate that at ambient 
conditions strong modifications of the water structure are 
induced by the solutes with the formation of hydration 
cages around the gas atom. Guillot and GuissaniO found 
that at increasing temperature the solute tends to dis- 
rupt the cage. The entropic contribution is dominated 
by the reorientation of the water molecules in the first 
shell of hydration, where only the repulsive part of the 



water-solute interaction plays a role. 

In this respect the change in the structural properties 
of water approaching the supercritical states with the 
modification of the hydrogen bond networktllE3E3 could 
be very relevant in determining the strong increase of the 
solubility of apolar gases with temperature. Recent neu- 
tron diffraction studies performed at supercritical condi- 
tions in mixtures of rare gases at low concentration in 
water found that the structure of water is modified and 
compressed by the presence of the apolar solutes.li-3 

Simple pair potential used in computer simulation of 
water based on parameters fitted to the properties at 
ambient conditions like TIP4PB and SPC/B3 do not 
agree with the experimental results of neutron diffrac- 
tion in predicting the ^urloiral changes in water at in- 
creasing temperatures-EjOlla Due to the interest in su- 
percritical water several polarizable potentials have been 
recently developed with the aim of reproducing th e pro p-, 
erties of water around and above the critical point oEj-EJ 
According to a recent detailed analysiso performed by 
comparing the p raalat ion with the more recent neutron 
diffraction dataEj-ESEa the polarizable models give pair 
correlation functions which agree with the experimen- 
tal data in the supercritical region much better than the 
nonpolarizable models. On the other hand the nonpo- 
larizable models reproduce better the temperature de- 
pendence of some thermodynamical quantitiesEj and the 
SPC/E potential for instance predicts a critical point for 
water in agreement with the experimentsEi, while some 
polarizable models do not reproduce well the coexistence 
curve £3 Polarizable models however are generally more 
accurate in describing the modifications of the structure 
of water at increasing temperature and pressure above 
the critical point .E3o 

In this paper we present a comparison of the pair cor- 
relation functions of mixtures of argon and water at su- 
percritical conditions obtained from computer simulation 
with the recent experimental data quoted above." In 
the numerical calculations water is simulated with the 
SPC/E potential and the polarizable potential model 
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proposed originally by Ruocco and Sampoli@ and further 
developed by Brodholt, Sampoli and Vallauri (BSV).Ej 
This model is based on the TIP4P geometry with a polar- 
izablc point dipole placed in the center of mass. This po- 
tential has been checked in computer simulation of pure 
water and found to produce pair correlation functions 
which are in good agreement with the experiment at su- 
percritical conditionsrjH In the next section we give 
some detail of our simulation. In Sec. 3 we present the 
results and compare them with recent experimental data. 
Sec. 4 is devoted to conclusions. 



TABLE I: Thermodynamical points where the computer sim- 
ulations have been performed. The point at T — 673 Kr-aor- 
responds to temperature and density of the experiment.tfl 



T(K) 
300 
373 
473 
673 



P (g/cm 3 ) 
0.998 
0.949 
0.841 
0.499 



II. COMPUTER SIMULATION OF AQUEOUS 
SOLUTIONS 

We performed the computer simulations of a mixture 
of 256 water molecules and rare gas atoms at the con- 
centration of one atom for 40 molecules of solvent, the 
same used in the neutron diffraction experiment quoted 
above. The interaction of the water molecules has been 
described by the polarizable BSV modeO and, for com- 
parison, by the SPC/E model.0 In both the cases the sim- 
ulations have been performed in the microcanonical en- 
semble with the minimum image convention and a cut-off 
of the interactions at half of the box length. The reaction 
held has been used to take into account the long range 
part of the electrostatic interactions. The system with 
the SPC/E water has been simulated with_the DLPoly 
package in the version 2.10, see referencetil for further 
details. 

The BSV potential has the same geometry as the 
TIP4P potentiate with a polarizable dipole moment 
placed on the center of mass of the molecule. The two 
positive charges, the parameters of the Lennard-Jones 
potential on the oxygen and the position of the negative 
charge are adjusted to reproduce some of the properties 
of water at ambient conditions (see model 4 in refo). 
The induced dipole pi = aEi is calculated from the local 
electric field E{ with an iterative procedure by assum- 
ing an isotropic polarizahility fixed to the value for water 
molecules a — 1.44 A 3 The long range part of the 
electrostatic interactions were taken into account by the 
reaction field method.E] Details of the extension of the 
reaction field method .to. include polarization effects can 
be found in referencesE3E3. 

We used a time step of 10~ 15 s and the structural quan- 
tities have been calculated by averaging over 2 • 10 5 con- 
figurations. We carried out the simulations for mixtures 
of water with Ar and Ne. Here we focus mostly on the 
results of the simulation of water with argon atops, for 
which we can compare with recent experiments.E-3 The 
interaction of the water with the noble gas atoms has 
been modeled with a Lennard-Jones potential between 
the oxygen of the water molecules and the solute. The 
parameters of the Lennard-Jones potential of argon are 
assumed as tAr jk B = 125 K, a Ar = 3.415 A and for neon 
^Ne/ks — 18.56 K and — 3.035 A.til The parameters 



of the oxygen-solute potential have been derived from the 
Lorentz-Berthelot rules. 

We explored a range of temperature and density from 
ambient conditions to the values in the supercritical re- 
gion where the experiment has been performed, as indi- 
cated in Table H. The calculations were performed below 
T = 673 at the points which correspond to the liquid, 
branch of the coexistence curve of the SPC/E model.E£l 
The point at T = 673 K corresponds to the conditions of 
temperature and density (« 0.499 g/cm 3 ) of the experi- 
ment. 



III. SOLUTE-SOLVENT STRUCTURE 

In Fig. [l]we report the pair correlation functions argon- 
oxygen goAr(r) at the different thermodynamical points 
indicated in Table | obtained by the BSV and the SPC/E 
models. At ambient conditions in both the models the 
goAr(r) show well defined first peak slightly different in 
the position due to the different value of the Lennard- 
Jones parameter. The argon-hydrogen pair correlation 
functions gHAr(r) (Fig. |^) show a first peak at the same 
position as goAr( r ) indicating that the solute is located 
interstitially in the hydrogen bond network equidistant 
on the average from the oxygens and hydrogens. 

The shoulder around 3.8 A in the gHAr(r) indicates 
however that some of the O-H bonds point radially to- 
ward the second shell. Comparing the positions of the 
second peak of the goAr{r) with the one of the gHAr(f) h 
is evident that the in the second shell the water molecules 
are preferentially oriented with the oxygen atoms towards 
the solute and the hydrogen atoms towards the bulk. 
Similar trends were already found bi£. Guillot and Guis- 
sani on rare gases in SPC/E watem and explained in 
terms of the formation of an hydration cage around the 
solute (seeB for more details). We can conclude that the 
same effect is predicted by the polarizable BSV model at 
ambient conditions. 

On increasing temperature the first and second peak 
of the goArir) decrease and the second hydration shell 
tends to disappear. In the gHAr{r) only the shoulder 
remains at the supercritical conditions indicating that 
the solute tends to expel the water molecules which are 
in its solvation shell. As it appears from Fig. [j] and || the 
BSV model confirms the prediction of the SPC/E model 
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FIG. 1: Oxygen-argon pair correlation functions on increasing 
temperatures: comparison between the BSV model (on the 
left) and the SPC/E model (on the right); temperatures from 
above are T = 300 K, T = 373 K, T = 473 K, T = 673 K. 
The thermodynamical conditions are indicated in Table EL 
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FIG. 2: Hydrogen-argon pair correlation functions on increas- 
ing temperatures: comparison between the BSV model (on 
the left) and the SPC/E model (on the right); temperatures 
from above are T = 300 K, T = 373 K, T = 473 K, T = 673 K. 
The thermodynamical conditions are indicated in Table f| 



and the microscopic scenario described in ref.EI, but it 
shows a more drastic reduction of the peaks in going from 
ambient temperature to the point at T = 373 K. 

In the supercritical state at T = 673 K the BSV poten- 
tial, which is more attractive with respect to the SPC/E 
model, predicts a more structured goAr{f) ■ We can fri- 
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FIG. 3: Oxygen-argon pair correlation function: comparison 
between the BSV model (solid line), the SPC/E model (long 
dashed line) and the experimental resultscj (crosses) at the 
thermodynamical conditions of the experiment (T = 673 K) 
as indicated in Table |l| 

fer that in the process of expelling the solvent the depth 
of the solute-solvent potential well plays an important 
role in agreement with the conclusions of conmitter sim- 
ulations!] and previous theoretical predictionsao on the 
relevance of the solute-solvent interaction as compared 
with the solvent-solvent. 

In Fig. the argon-water radial distribution func- 
tions obtained from oup-simulation are compared with 
the experimental resultaij. In both goAr{r) and gHAr( r ) 
the use of BSV model improves substantially the agree- 
ment with the experiment. The solute-solvent effective 
interaction results to be more attractive in the polariz- 
able model with respect to SPC/E as evident from the 
comparison of the height of the first peak of the radial 
distribution functions. 

This trend is confirmed by calculations of the coor- 
dination numbers of the water molecules around the 
argon atoms. The deficit of water molecules around 
the solute can be estimated by calculating the quantity 
h w (r) = [n Ar{r)/n u {r)\ - 1 where 

noAr(r) = 47175 / g Ar(r')r' 2 dr' (1) 
Jo 

and the coordination number for a homogeneous fluid is 
given by 

n u (r) = ^np(r 3 -a 3 ) (2) 

p is the water density and a is the minimum approach dis- 
tance between argon and oxygen. The results are shown 
in Fig. compared with the quantity h w (r) obtained from 
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FIG. 4: Hydrogen-argon pair correlation function: compari- 
son between the BSV model ( solid line), the SPrG/E model 
(long dashed line) and the experimental resultscfl (crosses) 
at the thermodynamical conditions of the experiment (T = 
673 K) as indicated in Table |. 
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FIG. 5: Deficit of water molecules around the solute, as de- 
fined in the text (see Eq. [l]), obtained in the BSV (solid line) 
and in the SPC/E (long dashed line) models compared with 
the same quantity calculated from the experimental points 
(crosses). 



FIG. 6: Radial distribution functions of the mixtures at the 
supercritical conditions of the experiment at T — 673 K (see 
Table |): gooir) solid line, goAr{r) long dashed line, gHAr{r) 
point line. On the left results frem the BSV model, on the 
right results from the experimented, in the inset results from 
the SPC/E model. 
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FIG. 7: Radial distribution functions of the mixtures from the 
BSV model at T = 673 K: goo(r) (long dashed line), goAr(r) 
(solid line), goNe(r) (point line). goo(r) is extracted from 
the Ar-water simulation and it is shown as representative of 
the position of the first O-O shell. 



the experiments. The negative values indicate a repul- 
sive behaviour of the solute, but the repulsion is overes- 
timated by the SPC/E model particularly in the region 
of the first peak of the goAr(r). 

In Fig. ^ we report together the radial distribution 
functions goo{ r ), 9oAr(r) and gHAr(r) for the water- 
argon mixture. On the right side there are the results of 
the BSV, on the left the experimental findings and the 



inset shows the SPC/E results. Both the models repro- 
duce correctly how the solute penetrates the shells of the 
water molecules. The goAr(r) does not penetrate beyond 
the second shell, while the minimum approach distance in 
the guAr (f) is shorter than that observed in the goo (j) ■ 
In Fig. ^ we show the solute-oxygen radial distribution 
functions of the water-argon and water-neon mixtures 
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FIG. 8: Radial distribution functions of the mixtures from the 
BSV model at T = 673 K: goo{r) (long dashed line), gHAr{r) 
(solid line), guNeir) (point line). goo( r ) is extracted from the 
Ar-water simulation and it is shown as representative of the 
position of the first O-O shell. 



together with the goo ( r ) of the argon- water mixture as 
obtained from the BSV model. The calculations for the 
water-neon mixture were performed at the same density 
and concentration of the water-argon mixture, which cor- 
respond, as stated above, to the experimental conditions. 
The goo ( r ) °f the two mixtures are very similar and any 
of them could be chosen for the comparison. It is evident 
the role of the size of the apolar solute. The Ne atoms 
almost penetrate the first shell of water molecules. The 
distance of minimum approach of Ne shifts to a lower 
value and roughly coincides with the minimum approach 
distance of the goo{ r )- The maximum of the goNe{ r ) is 
at the same position of the maximum of the goo{ r )- hi 
Fig. U it is shown that the gasi^) of each solute have a 
minimum approach shorter than that of the gos {r) ■ Our 



trend with the size of the solute is in agreement with the 
one found in the experiment. 



IV. CONCLUSIONS 

We performed a computer simulation of a solution of 
argon atoms in water at the same-conditions of a re- 
cent neutron diffraction experiments. We compared the 
results obtained by using .-the SPC/E modell and the 
BSV polarizable potential.Ej For both the models we de- 
termined the solute-solvent radial distribution functions 
from ambient to supercritical conditions. The local struc- 
ture of the solvent is perturbed by the presence of the 
solute with the formation of a cage of hydration around 
the non polar atoms. At increasing temperatures the 
solute tends to expel the water which resides in the sol- 
vation shell. Both the models reproduce this trend, at 
supercritical conditions however the BSV potential pre- 
dicts a more structured goAr( r ) with respect to SPC/E 
as a consequence of a more attractive solute-solvent in- 
teraction. The argon-water radial distribution functions 
are well reproduced at the experimental conditions by the 
polarizable model. The SPC/E model overestimating the 
solute-solvent repulsion gives the effect of collapsing the 
argon-water structure. The polarizable model predicts 
much better the solvent-solute correlation at supercrit- 
ical conditions. This supports the idea that-,thc BSV 
model works quite well at high temperatures .c3 
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